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FOREWORD 


This  report  was  prepared  by  Lt  T.  M.  Weeks  of  the  Electrogasdynamics  Branch,  Air  Force 
Flight  Dynamics  Laboratory,  Wright- Patterson  Air  Force  Base,  Ohio.  The  work  was  conducted 
under  Project  No.  1426,  “Experimental  Simulation  for  Flight  Mechanics,**  Task  No.  142604, 
"Theory  of  Dynamic  Simulation  of  Flight  Environment  (Dynamic  Flight  Simulation).* * 

The  period  of  work  covered  by  this  report  was  from  January  1966  to  June  1967.  This  report 
was  submitted  by  the  author  in  November  1967  for  publication. 


This  technical  report  has  been  reviewed  and  is  approved. 


PHILS  p.  ANTONATOS 
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ABSTRACT 


The  much  ignored  problem  of  turbulence  in  hypersonic  ground  test  f^oillties  is  discussed  with 
particular  emphasis  on  arc-heated  facilities.  Sources  of  free-stream  turbulence  and  their 
effects  on  model  aerodynamics  are  identified.  An  analytical  treatment  of  stagnation  cone 
heating  accounting  for  free-stream  vorticity  amplification  aud  viscous  dissipation  is  presented. 
Results  are  discussed  in  relation  to  three  experiments :  one  involving  time  resolved  measure¬ 
ments  in  an  arc-heated  wind  tunnel,  another  involves  stagnation  point  heat  transfer  measure¬ 
ments  in  the  turbulent  boundary  layer  of  a  hypersonic  shock  tunnel,  and  the  third  deals  with 
similar  measurements  in  a  turbulent  subsonic  nitrogen  plasma  jet.  hi  both  of  the  latter 
experiments,  stagnation  point  heating  rates  exceeding  two  to  three  times  the  calculated 
laminar  expectation  are  reported.  These  results  are  corroborated  by  the  present  analysis. 
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SECTION  I 
INTRODUCTION 


The  severe  reentry  environment,  both  present  and  future,  places  a  heavy  demand  on  ground 
test  facilities  to  provide  realistic  simulation.  A  problem  common  to  most  hypersonic  facilities 
which,  until  now,  has  been  virtually  ignored  concerns  the  presence  of  fre e-stream  turbulence. 
Broadly  defined,  hypersonic  wind  tunnel  turbulence  includes  all  unsteady  flow  behavior 
eventually  experienced  by  the  test  article.  Included  therefore  are  gross  property  variations 
in  the  test  affluent.  To  a  first  approximation,  frozen  disturbance  patterns  convected  by  the 
mean  flow  make  up  a  significant  portion  of  hypersonic  turbulence.  Inherent  in  the  nature  of 
hypersonic  turbulence  defined  in  this  way  is  the  lar  change  in  scale  of  a  given  disturbance 
in  passing  from  the  hypersonic  free  stream  to  the  model  flow  field.  Thus  free-stream 
disturbance  wavelengths  of  the  order  of  the  tunnel  size  may  become  of  classical  size  near  a 
stagnation  point. 

There  is  growing  evidence  that  such  hypersonic  turbulence  significantly  increases  both 
stagnation  region  heating  and  ablation  rates  as  well  as  reduces  the  transition  Reynolds  number 
when  comparing  these  effects  to  laminar  free-stream  results. 

This  report  discusses  the  general  aspects  of  the  overall  problem.  The  turbulence  is  charac¬ 
terized  and  an  anclvsis  of  stagnation  line  heat  transfer  based  on  earlier  work  of  Sutera 
(Reference  1)  is  presented.  The  effects  of  both  vorticity  amplification  and  viscous  dissipation 
arc  included.  An  Incompressible  flow  of  a  perfect  gas  is  assumed.  Some  experimental  results 
are  discussed  which  support  the  theory  and  provide  specific  evidence  of  the  importance  of  the 
problem. 


1 


AFFDL-TR-67-195 


SECTION  n 

THE  HYPERSONIC  TURBULENCE  PROBLEM 


1.  STATEMENT 

A  definition  of  hypersonic  wind  tunnel  turbulence  Is  given  in  Section  L  Any  spatiotemporal 
variation  of  any  flow  property  about  its  mean  value  is  a  component  of  hypersonic  turbulence. 
In  a  typical  ground  test  facility,  there  exist  certain  cause  and  effect  aspects  which  may  be 
treated  separately.  The  origin  of  the  turbulence  in  a  typical  hypersonic  facility  (excluding  the 
hypervelocity  range  facility)  is  clearly  related  to  the  particular  driver.  Figure  1  shows  an 
arc-heated  wind  tunnel  capable  of  providing  hypersonic  aerodynamic  simulation  or  low  Mach 
number  simulation  of  high  shear,  high  heat  transfer.  Several  mechanisms  which  may  render 
the  flow  turbulent  are  identified. 

Fluctuations  in  arc  heater  supply  voltage,  current,  and  pressure  are  obvious  sources  of 
turbulence.  If  the  arc  is  vortex  stabilised,  the  rotational  frequency  is  important.  The  heater 
Itself  may  be,  tor  a  certain  set  of  operating  conditions,  a  source  of  natural  longitudinal, 
radial,  and  circumferential  acoustic  resonance.  The  plasma  turbulence  in  the  cylindrical  cavity 
consists  conceivably  of  both  solenoidal  and  ir rotational  modes  as  described  by  Kovasznay 
(Reference  2).  The  task  of  identifying  various  sources  of  unsteady  flow  behavior  within  the 
heater  appears  formidable  owing  to  strong  coupling  effects. 

As  the  flow  is  accelerated  toward  the  throat,  additional  coupling  of  turbulence  modes  takes 
place  due  to  large  axial  and  radial  property  gradients.  The  fundamental  character  of  the  turbu¬ 
lence  changes  as  the  flow  becomes  supersonic.  The  various  separate  turbulence  modes 
examined  by  Kovasznay  (Reference  3)  and  Chu  and  Kovasznay  (Reference  4)  are  the  acoustic, 
vorticity,  and  entropy  modes,  ki  addition,  there  is  a  chemical  mode  which  in  arc-heated 
devloes  would  include  some  contamination.  Although  nonlinear  mode  cotpllng  still  occurs  near 
the  throat  and  within  the  boundary  layer,  one  expects  that  for  small  enough  disturbance 
amplitudes  the  various  modes  are  uncoupled  in  the  free  stream  (Reference  3).  If  the  nozzle 
boundary  layer  becomes  turbulent,  it  is  capable  of  producing  additional  free-stream  turbulence 
in  the  form  of  acoustic  radiation,  which  may  be  enhanced  by  probable  vibration  of  water- 
cooled  nozzle  walls. 

The  origin  and  basic  nature  of  the  turbulence  which  ultimately  appears  in  the  test  effluent 
have  been  characterized  above.  Any  one  of  the  described  property  variations  in  the  heater  can, 
by  the  mechanism  of  nonlinear  mode  coupling  near  the  throat,  give  rise  to  all  four  turbulence 
modes  in  the  test  gas.  For  this  reason,  all  hypersonic  facilities  operating  in  the  range  where 
chemical  kinetics  are  important  and  having  any  kind  of  disturbance  present  in  the  driver  are 
potentially  capable  of  containing  all  forms  of  turbulence  in  the  test  core.  This  first  aspect  of 
the  turbulence  problem  has  been  well  documented  for  supersonic  facilities  (see  Morkovin, 
Reference  5). 

The  second  major  aspect  of  the  problem  is  concerned  with  the  effects  of  the  various  free- 
stream  turbulence  modes  on  simulation.  The  three  separate  effects  identified  are  (1)  observed 
Increased  stagnation  region  heat  transfer  and  skin  friction  (References  6  through  12),  (2)  in¬ 
correct  mechanical  ablation  (Reference  13),  and  (3)  early  boundary  layer  transition  (Refer¬ 
ence  14).  The  ways  in  which  the  presence  of  free-stream  turbulence  brings  about  these  effects 
are  complex.  Basically,  each  of  the  turbulence  modes  (acoustic,  vorticity,  entropy,  chomical) 
is  capable,  qpon  shock  interaction,  of  generating  all  four  modes  behind  the  shock  (e.g.. 
References  3  and  14  through  19).  As  an  example,  fluctuations  in  free-stream  atomic  oxygen 
concentration  may  yield  fluctuations  in  velocity  behind  the  bow  shock  wave.  Behind  the  bow 
shock  we  may  again,  for  high  enough  stagnation  point  Reynolds  numbers,  identify  a  flow  region 
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ahead  of  the  boundary  layer  herein  the  four  modes  are  uncoupled.  Convection  of  vorticity 
enthalpy  and  chemical  modes  and  propagation  of  acoustic  waves  through  the  stagnation  point 
boundary  layer  will  alter  the  surface  heat  transfer  rate  and  skin  friction  calculated  by 
assuming  laminar  flow  and  laminar  transport  properties  (including  viscosity,  conductivity, 
and  species  diffusion  coefficients).  The  role  of  wall  catalycity  must  be  reaxamined,  and 
surface  catalytic  efficiency  will  be  seen  to  depend  on  boundary  layer  turbulence  levels.  Tur¬ 
bulent  flow  in  the  boundary  layer  will  also  quite  naturally  result  in  modified  ablation  patterns. 
Early  transition  to  a  turbulent  boundary  layer  on  the  model  bears  an  obvious  connection  to  the 
free-stream  turbulence  level  found  in  conventional  supersonic  wind  tunnels. 

It  is  not  Intended  to  suggest  that  each  of  the  three  effects  exists  separately  because  all 
three  might  be  encountered  in  a  typical  hypersonic  wind  tunnel.  The  major  point  is,  however, 
that  for  both  cited  applications  of  the  arc-heated  facility,  i.e.,  hypersonic  aerodynamic  and 
high  shear,  high  heat  transfer  simulation,  the  problem  of  free-stream  turbulence  must  not  be 
overlooked. 

2.  APPROACH 

The  first  of  two  possible  independent  approaches  to  the  problem  of  hypersonic  wind  tunnel 
turbulence  is  a  direct  attack  aimed  at  eliminating,  suppressing,  or  at  least  altering  the 
character  of  the  turbulence  sources.  The  usual  turbulence  suppression  devices,  well  under¬ 
stood  in  supersonic  wind  tunnels,  cannot  be  installed  in  hyperthermal  facilities  without  serious 
loss  of  total  enthalpy.  Of  course,  a  high  degree  of  rectification  of  the  power  supply  and 
“compressor  off*  operation  are  indicated.  Careful  tuning  of  the  heater  may  reduce  acoustic 
resonance  effects  for  a  particular  operating  point.  Fluctuations  in  arc  length  and  attachment 
mode,  a  major  contributor  to  effluent  turbulence  especially  at  high  driver  pressures,  are 
still  without  apparent  solution. 

Another  important  suppression  technique  might  be  aimed  at  reducing  property  gradients  in 
the  throat  region.  The  present  day  cone  cyclinder-cone  throat  geometry  would  probably 
require  modification.  A  suggestion  based  on  experience  with  low  power  level  choked  flows 
indicates  that  uncooled  graphite  throats  yield  much  quieter  flows  than  highly  cooled  copper 
throats  (Reference  20).  It  would  appear  that  the  greatly  reduced  radial  temperature  gradient 
in  the  graphite  throat  case  reduced  noise  generation  by  the  coupling  mechanism  previously 
described  in  this  Section. 

hi  order  to  bring  about  a  reduction  of  arc-heated  wind  tunnel  turbulence,  it  will  be  necessary 
to  initially  characterize  the  various  modes  of  turbulence  which  are  present.  This  may,  in 
part,  be  accomplished  via  a  number  of  flow  diagnostic  techniques  currently  available,  although 
by  no  means  well  understood  in  the  hyperthermal  range.  Thus  one  is  led  to  consider  the  various 
flow  visualization  techniques  including  laser  doppler  schlleren  and  crossed  optical  beams 
(Reference  21),  as  well  as  inserted  probe  devices  such  as  hot  wire  and  cooled  sensors 
(References  22  and  23),  and  acoustic  and  heat  flux  probes  (see  Section  IV).  A  promising  tech¬ 
nique  employing  an  electron  beam  with  fiber  optics  has  been  used  successfully  (Reference  24). 

The  second  approach,  which  can  be  carried  out  in  large  part  independently  of  the  first, 
involves  the  analysis  of  the  effects  of  arbitrary  forms  of  turbulence,  characterized  for  example 
by  amplitude,  frequency,  and  correlation  on  model  aerodynamics.  Major  emphasis  should  be 
given  to  whatever  modes  of  turbulence  are  actually  present  in  an  arc-heated  facility.  At  the 
present  time,  it  would  appear  safe  to  assume,  based  on  previous  arguments,  that  all  modes, 
albeit  uncoupled  in  the  fires  stream,  are  present.  The  presence  of  these  free-stream  turbu¬ 
lence  modes  results  in  modified  stagnation  zone  flow  which  may  be  subdivided  into  two 
categories:  (1)  shock-turbulence  interaction  where  the  newly  created  or  modified  turbulence 
modes  behind  the  bow  shock  must  be  identified  and  (2)  the  model  boundary  layer  flow  field 
having  initially  uncoupled  turbulence  modes  at  its  outer  edge.  The  order  of  attack  is  not 
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important,  and  each  major  component  maybe  subdivided  into  various  mode  coupling  problems, 
such  as  shock- vorticity  (References  5  through  8);  shock-chemical  mode  coupling;  boundary 
layer-vorticity;  boundary  layer-acoustic  mode  coupling. 

Recently,  Sutera  (References  1  and  25)  analyzed  the  effects  of  free-stream  vorticity  ampli¬ 
fication  through  a  stagnation  line  boundary  layer  on  skin  friction  and  heat  transfer.  He  found 
that  both  the  skin  friction  and  heat  transfer  increased  (by  different  amounts)  due  to  the  pres¬ 
ence  of  a  specifically  oriented  upstream  steady  vorticity  field  convected  through  the  boundary 
layer.  The  analysis  was  carried  out  assuming  incompressible  flow  of  a  perfect  gas.  Provision 
was  made  to  include  the  generation  of  higher  harmonic  spectral  interaction  functions  including 
temperature- velocity  coupling  terms.  Effects  of  viscous  dissipation  were  not  considered 
and  the  actual  numerical  solution  was  accomplished  using  truncated  forms  of  the  higher  order 
approximations  because  of  the  limited  analog  computer  capacity. 

The  effects  of  species  diffusion  should  be  included  although,  as  stated  previously  in  this 
Section,  this  analysis  should  be  based  on  turbulent  diffusion  coefficients  as  yet  not  well  known. 
R  would  appear  then  that  a  logical  extension  to  Sutera’ s  problem  is  to  assess  the  effects  of 
viscous  dissipation  and  include  all  higher  harmonic  terms  in  the  highest  practical  approxi¬ 
mation.  This  represents,  approximately,  the  flow  field  behind  a  hypersonic  bow  shock.  One 
may  regard  the  presence  of  the  spatially  distributed  vorticity  field  as  a  direct  result  of  shock- 
turbulence  interaction. 
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SECTION  III 

TWO-DIMENSIONAL  STAGNATION  POINT  HEAT  TRANSFER 
WITH  VORTICITY  AMPLIFICATION  AND  VISCOUS  DISSIPATION 


1.  MATHEMATICAL  STATEMENT  OF  PROBLEM 

The  need  to  include  effects  of  viscous  dissipation  arises  from  the  result,  obtained  by  Sutera 
et  al  (Reference  25),  that  a  neutral  scale  vorticity  exists.  For  turbulence  wave  lengths  above 
t'lJLs  scale  the  mechanism  of  vorticity  amplification  is  found.  For  wavelengths  below  this  scale 
there  is  vorticity  attenuation.  Sutera  attributes  this  to  viscous  dissipation  but  does  not 
formally  provide  for  this  effect  in  the  energy  equation. 


In  this  report,  the  basic  approach  to  the  problem  initially  parallels  that  of  Sutera,  i.e.,  the 
same  basic  form  of  the  upstream  vorticity  with  wavelength  again  equal  to  1.5  times  the  neutral 
is  chosen  for  purposes  of  comparison.  Figure  2  shows  the  basic  flow  model.  The  added 
vorticity  is  in  the  only  direction  capable  of  amplification.  The  modified  analysis  begins  with 
Sutera’ s  form  of  the  energy  equation  but  with  the  viscous  dissipation  term  included. 
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Specializing  to  the  Hiemenz  flow  configuration  which  employs  nondimensional  flow  quantities 
and  the  same  scale  factor  one  has: 


V  )  T  = 


Pr 


V2  T 


QV 


cp  -Tw> 


4> 


(2) 


The  quantity  Jov  has  been  used  to  nondimens ionalize  the  velocity  and  one  recognizes  the 
dimensionless  Eckert  number: 


E  = 


QV 

Cp  (T^  —  Tw  ) 


(3) 


Note  also  that  the  quantity  v/a  has  the  dimensions  of  length.  The  complete  system  of 
governing  differential  equations  expanded  into  the  £ ,  rj ,  £  nondimensional  coordinate  system 


becomes  (Reference  1): 

-&U  =£1^1)  +  &££  (4) 

UUC  +  Vuiv  +  WUCC  “u5%  +  uv  VC  *  “V?  +  UK  >5  15) 

uu,  +vu^,  +wu^  +  uc  w,-  u^wc  =  (Uvn  +u5C  )v  ,6) 

U  +  V,,  +  Wf  =0  (7) 

VTv  +  WTc  =  ■? r  ,Tw  +Tcc '  + E  •  *  (si 

where 

4>  =2  (  u2  +  v2  +  w2  +  )  +  w2  -P  v^2  (9) 


2 

S  has  been  assumed  in  Equations  4  through  9  that  terms  involving  £  in  the  dissipation 
function  could  be  neglected.  This,  of  course,  restricts  the  solution  to  a  region  close  to  the 
stagnation  line.  Next,  it  is  assumed  that  the  three-dimensional  nature  of  the  turbulent  flow 
field  may  be  characterized  by  the  following  Fourier  structure  which  permits  solutions 
periodic  in  £ . 

U  =  U0  (77)  +  A  £  «n  co*  kn  5 
n  =  l 

• 

V  =  V0  ( r 7)  +  A  £  vn  (77) cos  kn  £ 

.  n=l  (10) 

W  =  I  -j-  k  £ 

n  =  l 

T  =  T0  {7j)  +  A  I  6n  (77)  cos  kn  £ 
n  =  l 
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The  harmonic  coefficients  uQ,  vq,  q  q  are  not  necessarily  small  with  respect  to  their  mean 

flow  counterparts  UQ,  VQ,  T^,  respectively.  Substitution  of  Equations  10  into  Equations  4 

through  8  leads  ultimately  to  a  set  of  two  plus  a  quintuple  infinity  of  ordinary  nonlinear 
differential  equations  with  as  many  unknowns.  These  equations  express  both  mean  and 
"fluctuating*  *  quantities  as  functions  of  distance  from  the  wall.  The  detailed  set  of  equations 
without  viscous  dissipation  appears  as  Equations  18  through  27  of  Reference  1.  The  energy 
equation  governing  the  mean  temperature  variation  modified  to  include  dissipation  becomes: 

T0"  +  *  <f>  T0*  +  4E  Pr  <f>a  =  Pr  I  [Uj  0,  +(vj  0j  )'  ] 

*"1  (II) 

-A2  Pr  E  I  [2  (u*  +w2  +v!2)  -2  v.  w.'  +  (v.  k,)24-  (Wj'/k.  )  ] 
i*  I  1 

The  fluctuating  temperature  field  is  governed  by  the  following  egression: 


»n+P'  *  S:-9n  *1  '  P'  'o'  »,  +4E  Pr  f  (u„  -v„'  )  •  f-  Pr  » 


-AE  Pr  { I_  +  I  +  n  +  I-  _+  I.  n) 

3,n  4,n  5/1  6,n  7,n  s,n 


(12) 


The  Spectral  Interaction  functions  I  n  ^  are  tabulated  in  the  Appendix.  It  may  be  seen  that 

even  without  free-stream  turbulence  (A=o)  there  is  a  viscous  dis station  effect  owing  to  the 
finite  value  of  UQ  *  -V^  *  <£.  Solutions  of  the  entire  system  of  equations  have  been  obtained 

by  means  of  an  analog  computer  and  a  digital  computer  simulation  of  the  analog  computer  pro¬ 
gram. 

2.  ANALOG  COMPUTER  SOLUTION 


S  is  necessary  to  choose  a  wavenumber,  kj ,  for  the  first  harmonic  component  of  the  initial 

velocity  and  temperature  field.  Sutera  has  considered  a  simple  sinusoidal  variation  of  the 
normal  velocity  which  corresponded  to  a  single  value  of  k^  at  the  outer  edge  of  the  boundary 

layer.  One  can,  in  principle,  include  additional  wave  numbers  to  more  closely  represent  a 
hypothetical  turbulent  velocity  distribution.  Such  a  modification  would  not  alter  the  fundamental 
heat  transfer  mechanisms  under  examination  and  would  unnecessarily  complicate  the  analysis. 
The  appearance  of  higher  harmonics  within  the  boundary  layer  is  included  in  the  analysis  and 
it  is  seen  that  only  the  second  harmonic  terms  must  be  retrained. 

The  most  appropriate  wave  number  to  choose  is  one  which  exceeds  the  neutral  scale  value 
thus  ensuring  that  both  the  vorticity  amplification  and  viscous  dissipation  mechanisms  will 
be  effective.  It  is  necessary,  therefore,  that  k^<  1  which  is  the  same  as  X,  >  X0  where 

XQ  is  the  neutral  wavelength.  To  facilitate  direct  comparison  with  Sutera’s  result,  a  value  of 

kj  *  2/3  was  chosen.  A  range  of  the  turbulence  amplitude  parameter  A  was  considered  for 

several  values  of  the  stagnation  zone  Eckert  number.  The  latter  values  were  chosen  about  the 
usual  narrow  range  encountered  in  hypersonic  facilities. 

The  iterative  solution  began  by  assuming  that,  to  first  approximation,  the  spectral  interaction 
functions  contribute  nothing  to  the  velocity  and  vorticity  components.  The  homogeneous  forms 
of  the  governing  equations  were  solved  simultaneously  and  provided  a  first  approximation  to  the 
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first  and  second  harmonic  components  (n=l,2)  of  velocity  and  vorticity  interaction  terms. 
These  components  were  then  substituted  into  the  governing  equations  to  yield  a  second  approx¬ 
imation  to  the  velocity  and  vorticity  field.  This  procedure  'generates  a  second  approximation 
to  both  the  first  and  second  harmonic  spectral  interaction  functions.  These  results  were 
applied  to  the  energy  equations  (both  mean  and  harmonic  forms)  to  obtain  first  approximations 
to  the  spectral  temperature-velocity  interaction  functions  which  finally  yield  values  of  the 
mean  temperature  field. 

The  actual  calculation  scheme  as  it  was  programmed  for  the  analog  computer  and  sub¬ 
sequently  written  in  MIMIC  (digital  simulation)  is  presented  below.  The  set  of  homogeneous 
differential  equations  governing  the  first  approximation  to  the  first  harmonic  velocity  and 


vorticity  components  are  listed  first 

+M  *%)*{>)' +M+  1=0  "3) 

(See  Schlichting,  Reference  26) 

-V.'(.)  +'f'Vi(.>_  “f"  *■(.)“  u.'(tJ  =  0  04) 

ul'(i)  +<fyi)ui(i)  -“•(»)  (2<^(i)  ]  ~<f>U)  vi(i)  05) 

*.(.)•  "v.'(.)  ”  U.C.)  06) 


The  corresponding  boundary  conditions  to  be  Imposed  are: 

•*>(.)«»  =  ^.)(0)  =  °  ’  <*,(',){0)  =  l  237 
v.(.)(0)  =  V(.)(0)  *° 

ut(,,(0)  =  0 

4>\ ,)  ~ >  08  v  -* 00 
Mo  0  08  7?~>G0 

A  boundary  search  routine  was  employed  in  the  MIMIC  program  and  the  solution  to  Equa¬ 
tions  13  through  15  was  obtained.  Spectral  velocity  interaction  functions  were  defined: 

"N* 

^2(1 )  =  “Ml)  W|(.)  +  aMi)vi(i)+2w((i)<*M.) 


and 


U 


2  ( I ) 


uk.  > 1 
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Then  the  governing  equations  for  a  first  approximation  to  the  second  harmonic  (k^  =  4/3) 
velocity  (and  vorticity)  components  become: 


n  ,  16  8  '  - 

“v*(0  +  T  v2(i)  T  "eCi)  u2{  i)  =0 


W8(l)  +  <fyl)  U,2(l)4'^(l)  Wt(.)  "  9  W2{llS  2  ^2(1) 


«*|.1  +  <£m  u 


8  (l )  T(.)  U2(,)-U8(,)  ^(,)T9 


(|)V8(|)  '  2  W8(1) 


Uf  S  *"•  tf  ^  II 

*(  I  )  2(|)  U8(|) 


The  boundary  conditions  are: 

W0)  *  Vj(,|(0)  =  U2<,)«» 

•ill)  •  vf(.)  ■  U2(|>  ->0  °*  V  -*■  ® 

Solutions  of  this  set  of  equations  were  utilized  to  obtain  the  following  interaction  functions: 

u r “"<(■)  v8(.r w»(i)  v*(»r  wi(ow8(i)+a,8(»>  v!(ii+wawv.(i)"2*  w*o)  w'(«) 


Ul(8)=  4ul(.)U2(.  )+  Ml  )  U  2  ( I  )  +V./(I)U8(,)-'  W.(.  )  U8(.) 

V2  (i )  ui(l)+  V8(l)  Ul(l)+  ~Z~  W8(l)  Ul(l) 

The  next  step  was  to  calculate  the  second  approximation  to  the  first  harmonic  velocity  and 
vorticity  field. 

U I  ( 2)  ~  M2)  +  -f-  vl( 2) ~  “f-  ".(8)  =0  (20 

wi(8)  +  <fyl)U\(2)  +^(0W|(8)  "9"  Wl(a)  =  ^1(8)  (22) 


u 11  •+  u  *  —  u  ( 2 cb  +  ~  )  •  <f>*  ,  y 

1(8)  M)  1(2)  1(2  )’  ^  '  (l)  ~  9  '  ”(l)  V 

W.(8)  S  Vl(2)  ~  U. (8) 


fl)  V((2)=  2  ^1(2) 


The  boundary  conditions  are: 


'.<») (0>  1  Vui<°>  =  *  0 

*■'(«)•  ,.<2)  •  °S  V  -*a> 
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The  final  velocity-vorticity  interaction  functions  obtained  were: 

^2(2)=  "1(2)  V«'(2)+  a,l'(2)  V|(2)  +  2wi(2)a,|(2) 


U2(2)  =  2(u|(z)+  wl(2)  ul{2)  )  4  (v|(2)ul(2)) 

A  correction  function  for  4>  representing  contributions  of  the  spectral  interaction  function  to 
the  mean  velocity  field  was  also  calculated. 

%  =uf(2)+~2‘  { ui(z)v!(2}+ui(2)vi(2p+  4(.)+T(u*Ci)v*(0+u*(.)W 

which  yielded  a  second  approximation  for  <j> ,  i.e., 

<t>"  4  d>  d>‘  -  d>'Z  4  I  =  A2  ^  (25) 

■(2)  M2)  M2)  ^(2)  A  (2 ) 

The  final  system  of  equations  solved  for  the  velocity  and  temperature  field  corresponded  to 
the  second  approximation,  second  harmonic  components. 

'  «  .  J6  -o  (26) 

U2(2)“  v2(2)  +  9  v2(2)  3  W2{2)  ° 


a,2(2)4-^(2)U,2(2)+  4)<U2(2)"*9'  W2(2)=  f“  ^«<*> 


U2(2)  +^(2)U2(2)  U2(2)(2<^(2)+  9*  ^(2)  *2(2)  =  2  ^2(2) 


W2(2)=_W2(2)  U2  ( 2) 
The  boundary  conditions  are: 


V2(2)(°)  =v2(2)<°>  =u2(2 /0)  =0 
v2(2)i  v2(2)i  U2(2)— ►  0  0*  7)  — >  ® 


A  third  approximation  to  4*  was  obtained,  i.e.. 


<b  4  d>  —  fb'  4  I  =  A  ^ 

Ms)  Ms)Ms)  y(s) 


where 


^(3  )  =  Ul(2)4'  2  (U|(2)V|,{2)+  °l(2)  Vl(2),+  U2(£)4  2  (  U2(2)V2(2)  +  °2  (  2)V2(2)  * 
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The  system  of  Equations  11  and  12  governing  the  mean  and  fluctuating  temperature  required 
several  spectral  velocity  interaction  functions  defined  as  follows: 


®(,V=  8  [UI(2)U2  (2)+  V|(2)  V2  (2)4  W2(2)w|(2)]  +  9  V|(  2)  V2  <Z) 


+  4  W2{2)Wl{2)  W2  (2) Vl  (z)  ~  4  Wl  (2)V2(2) 


(31) 


®(2)=4[U"(2)+  V.U)  +*.%)] 


9  ,2 


9  l(2>  16 


*1(2)  +  2  Wl'(2)  Vl(2) 


(32' 


2  2  2  2  2  2 

T(  2 )  =  2(U,(2)+W,{2)+  V'(2)+U2(2)-hW2(2)4-V2(2)) 


.  I  ,  2  3  i  ,z  ,  I  .  4  3  /  ,2 

+  2  3  Vi  ( 2 )  2  w  I  (2)  +  2  3  V2(2)  4  W2(2)* 

The  first  approximation  to  the  mean  temperature  field  is  obtained  from  Equation  11. 

<2 


(33) 


with  boundary  conditions: 


Pf  T%)^3)+<(,)+4EPr  4»=0 


T0.  (0)  =0 


(34) 


'(■) 


'(.) 

;  'o(|).  T0"(|)  ->0  OS  7 7  -><D 


The  first  approximation  to  the  first  harmonic  component  of  the  fluctuating  temperature  field 
is  obtained  from  Equat.-on  12  which  becomes: 


®l'(l)+Pr  ®,(l>  •>>(!)-  Pr  T0(i),1(2)+4E  Pr  ♦'(*)  [  »f(*r  "l(«l 


A  Pr  • 

+  -5- (I,  -I,  >  -AE  Pr  I  I, 

2  (1,1)  (2,i)  Q^3  (or,  1 ) 


(35) 


The  spectral  interaction  functions  I  are  listed  in  the  Appendix,  Appropriate  boundary 
conditions  take  the  form;  '  ’  ' 

01(|)(°)  =0  ;  6lU)  ->  0  os  7j 


00 


The  first  approximation  to  the  second  harmonic  component  of  the  fluctuating  temperature 
field  may  be  evaluated  from: 


e,(„  +  pr  -  -f-  9?(l,=  Pr  To(1jV2(2,+  4E  Pr  [  v 2,- u2  (a, 


+  ¥r[I(,,zfI<d'AEPr  X  I(«' 


(36) 


2) 
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with  boundary  conditions: 

«„l)  <0)  =o;  (0  * 0  05  v  - 

We  are  now  in  a  position  to  evaluate  the  second  approximation  to  the  mean  temperature  field. 


oo 


TV, +  p'  t»I21  +  4E  Pr  *m  -  T*  [ 1 .  4  E XJ 


where 


'(2)  r(3)  "(2)  ''(3)  2 


*1  =  U|(2)^|(|  )  +  ^V|(2)^|(|p  +  U2(  2)^2  (I)  ■*"  ^  V2  (2)^2  (l)  ^ 


(37) 


V  2  <“l(2,+  ,;(2>  +,1(2)>-2,,(2)",'(2)+  T  "M  *  7  *.'(2) 

+  2<U2(2)  +  ,2(2)  +W2U)>-2V2(2)W2<04  V  *2(2)  +  l?  ",'(2 1 

Solution  of  Equation  37  provides  necessary  input  to  the  calculation  of  the  second  approximation 
to  the  first  harmonic  component  of  the  fluctuating  temperature  field  which  is  obtained  from  the 
relation: 

Sl'(2)+P'  ®,  (2)  *(3)-  -5-0.<2)  *P'  T0U1».(«]+''E  Pr  *(,).[«,' (*)-  “,(2)] 

,06) 

+  V(Vo-i(2,,)>-aep'  ^  I(a.») 

a=e 


The  second  and  final  approximation  to  the  second  harmonic  component  of  the  fluctuating 
temperature  field  may  be  obtained  from: 


K'(t)  +Pr  @2(2)  ^2(2)=  Pf  T0,,,V2(2)+4E  Pt  ( V2  (2)”  U 


'(2)  zi *■)  2(2 ) 

APrr  1  J£ 

+  ~5 —  I.  .-I.  .  -AE  Pr  I  I( 

2  L  (i,2)  (2,2)J  _  ta. 


(39) 


2  ) 


(40) 


The  results  obtained  above  may  be  combined  to  provide  a  thfrd  and  final  approximation  to  the 
mean  temperature  xield. 

T»(„+P'  Pr  <,  ■  ^T~  l1-  +  EI*1 

Where,  consistent  with  the  utilization  of  the  best  available  approximation,  we  have: 

II  =  Ul(2)^l(2)-f  '  vl(2  )  (2  )  >  +  UZ(2)S2(2)+  {v2(2)^2(2)} 

and  Ig  remains  unchanged. 
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3.  RESULTS  OF  ANALYSIS 

The  important  results  of  the  analysis  are.  of  course,  the  effects  of  vorticity  amplification 
by  vortex  stretching  and  of  viscous  dissipation  of  vorticity  on  stagnation  line  heat  transfer. 
S  is  worthwhile  in  this  connection  to  examine  the  intensity  amplification  of  the  first  and  second 
harmonic  components  of  v  and  u  referred  to  the  corresponding  local  mean  velocity.  Thus  we 
may  define  the  “turbulence  intensity”  ratios. 


II 

c 

r* 

/ 

U0 

=  1  Ul(2)  1 

/ 

*'<3, 

1,2  =  I  Vl  (  2  ) 

1  / 

V0 

=  1  *.  (,) 

1/ 

^(3) 

2,1  =  |  U2  (2) 

1/ 

*'<3> 

2,2  "1  W2  (2) 

1/ 

*(3) 

These  ratios,  plotted  in  Figure  3  as  a  function  of  y,  depict  the  way  in  which  the  harmonic 
components  of  the  disturbance  velocity  field  are  amplified  over  and  above  the  corresponding 
local  mean  velocity  as  the  flow  approaches  the  wall.  This  behavior  is  in  qualitative  agreement 
with  the  results  of  Reference  27.  Included  is  Sutera’s  (Reference  1)  result  for  r  *  J  /VQ. 

One  sees  that  Sutera  predicts  a  relative  intensity  amplification  which  exceeds  unity  and  also 
displays  a  maximum  for  V  >  O.  The  present  result  predicts  apeak  intensity  of  0.8  at  V  =  O. 
One  may  still  conclude,  however,  that  there  is  of  the  order  of  a  hundredfold  amplification  of 
initial  longitudinal  velocity  intensity.  The  first  harmonic  component  of  the  lateral  intensity 
attains  a  much  smaller  magnitude  (a  maximum  of  4%).  The  corresponding  second  harmonic 
intensities  are  also  of  smaller  magnitude  than  the  longitudinal  first  harmonic.  The  longitudinal 
second  harmonic  intensity  displays  a  double  maximum  owing  to  a  sign  change  in  ^2(2)' 

The  combined  effects  of  vorticity  amplification  and  viscous  dissipation  may  be  understood 
by  examining  the  source  term  in  Equation  40.  Thus  we  may  define: 

N*  =  [l,  +  EI,]-4E  P-  Ml) 

This  reduces  in  form  to  Sutera’s  for  A-l,  K=0.  Plots  of  N*  vs  V  appear  as  Figure  4.  In 

the  case  of  laminar  flow  (A=0)  we  find  that  the  source  term  is  non-zero  for  E  >  O.  One  may 
examine  the  variation  of  N*  as  a  function  of  A  alone  only  for  the  case  that  E=0.  For  E  >  O  one 
cannot  separate  the  two  effects  owing  to  the  middle  term  in  Equation  41. 

Solutions  of  Equation  40  in  the  form  of  T  ,  T '  are  plotted  vs  rj  in  Figure  5  for  the 

°<3)  °<3) 

same  range  of  A  and  E  as  appears  in  Figure  4.  The  overall  effect  of  increasing  E  is  to  in¬ 
crease  the  wall  temperature  gradient  and  reduce  the  thickness  of  the  thermal  boundary  layer. 
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Figure  4.  Source  Term  in  Energy  Equation  (Equation  40)  Inciuder  Both  Effects 
of  Viscous  Dissipation  and  Vorticity  Amplification 


Figure  5. 


Mean  Temperature  and  Temperature  Gradient  Illustrating  Distortion 
Due  To  Viscous  Dissipation  and  Vorticity  Amplification 
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The  level  of  T  labeled  “S”  Is  Sutera’s  value  of  T  .  This  level  which  is  to  be  compared 
(3)  °  (2) 

with  the  E=0  solution  is  somewhat  higher  in  all  cases,  the  difference  increasing  with  the  value 
of  A. 

We  come  finally  to  the  most  significant  result,  namely  the  wall  heat  transfer  rate,  q  which 
is  proportional  to  T*  (O).  For  a  reference  q.  we  chose  T  '  (O)  for  A  *  E  =  O  which  has 

°<3)  °(3) 

a  value  of  0.499.  Dividing  other  wall  values  of  T'  by  this  value,  we  can  show  in  Figure  6 

0  (3) 

that  the  effects  of  both  vorticity  amplification  and  viscous  dissipation  are  to  significantly 
increase  stagnation  line  heat  transfer.  As  previously  indicated  in  this  Section,  owing  to  the 
large  gradient  in  V  near  the  wall  there  is  finite  viscous  dissipation  even  for  laminar  flow. 

One  would  expect,  therefore,  that  even  in  free  flight  an  underestimation  of  stagnation  heat 
transfer  at  a  leading  edge  will  arise  (the  level  depending  on  the  magnitude  of  E)  if  viscous 
dissipation  is  not  accounted  for.  One  observes  that  the  effect  of  increasing  A  (as  noted  by 
Sutera)  or  E  is  to  increase  the  stagnation  zone  heat  transfer.  The  effect  of  viscous  dissipation 
is  initially  linear  but  is  gradually  reduced  as  the  boundary  layer  becomes  thinner.  This  is  due 
to  the  fact  that  the  source  terms  N*  In  Equation  40  (see  Figure  4)  while  themselves  extending 
well  beyond  the  classical  Hiemenz  boundary  layer  have  reduced  influence  on  the  thermal 
boundary  layer  which  grows  thinner  with  increasing  E.  This  effect  was  observed  by  Sutera 
for  Increasing  Prandtl  number. 

Included  in  Figure  6  are  four  calculated  values  for  the  Eckert  number  corresponding  to  the 
stagnation  flow  region  of  both  an  air  arc-heated  wind  tunnel  where  =  1.2  and  a  perfect  gas 

facility  where  Y  -  1.4.  fo  two  cases  (a,  b)  the  hypersonic  limiting  value  of  Mj  «  0.4  was 

chosen.  In  the  third  case  (c)  a  free-stream  Mach  number  =  2,  =  1.2  was  selected.  The 

final  case  represents  a  conventional  supersonic  wind  tunnel  condition  with  Y  -  1.4,  *  3. 

Transonic  and  subsonic  wind  tunnels  would  appear  with  still  higher  values  of  E.  In  all  cases 
the  Eckert  number  was  obtained  from  the  relation 

E  =  {  y  ~  I  5  M2) 

The  values  obtained  checked  well  with  Eckert  numbers  calculated  directly  from  the  definition 
(Equation  3). 
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SECTION  IV 

EXPERIMENTAL  RESULTS 


1.  SHOCK  TUNNEL  TURBULENT  NOZZLE  BOUNDARY  LAYER 

Some  recent  measurements  by  J.  E.  Wallace  (Reference  28)  of  stagnation  point  heat  transfer 
obtained  with  axisymmetric  probes  In  a  shock  tunnel  nozzle  wall  boundary  layer  provide  a 
comparison  with  the  theoretical  predictions  obtained  in  Section  IQ  S  was  found  that  “much 
higher  heating  rates  occurred  in  the  boundary  layer  than  would  be  precicted  based  on  measured 
pitot  pressure  profiles  and  an  assumed  Crocco  energy  relation  for  the  total  enthalpy  varia¬ 
tion.”  The  reason  was  attributed  to  the  tree-stream  turbulence  described  in  References  1,  6, 
and  7. 


fc  order  to  construct  such  a  comparison  it  is  necessary  to  relate  the  parameter  A  to  the 
initial  turbulence  amplitude  |  v  |  /V  ahead  of  the  probe  boundary  layer.  Since  the  two  are 

linearly  related  we  need  only  to  establish  a  single  data  point  where  both  the  excess  heating 
over  and  above  the  laminar  flow  value  as  well  as  the  free-stream  longitudinal  turbulence 
level  are  known.  Of  course  one  would  expect  that  a  meaningful  comparison  must  necessarily 
be  limited  to  two-dimensional  configurations.  Nevertheless,  one  may  anticipate  that  the  proper 
order  of  magnitude  effects  for  axisymmetric  shapes  would  be  characterized.  It  is  more  than 
likely,  however,  that  the  constant  of  proportionality  between  A  and  |  v|/V  is  different  for 
the  two  cases.  ° 

Short  et  al  (Reference  9)  reported  measurements  of  both  free-stream  turbulence  and 
corresponding  heat  transfer  rates  from  a  0.5  inch  diameter  sphere  in  a  low  speed  (16  fps) 
flow.  They  reported  that  an  overall  heat  transfer  increase  of  20%  accompanied  a  free-stream 
turbulence  Increase  from  0.Q13  to  0.15,  If  we  assume  that  at  the  lower  reported  turbulence 
level  the  heat  transfer  corresponds  nearly  to  the  laminar  rate,  we  find  from  Figure  6  a  value 
of  A  *  3.8  corresponding  to  E=0  and  q/q^  =  1.2.  This  provides  a  relation  between  A  and 

|  v|/VQ  for  axisymmetric  bodies  in  subsonic  flow  of  the  form: 

a=(t£>M/v 

For  the  shock  tunnel  experiments  it  is  again  necessary  to  know  the  turbulence  level  ahead  of 
the  probe  boundary  layer.  This  may  be  related  to  the  longitudinal  velocity  turbulence  level 
(axial  component)  actually  present  in  the  supersonic  boundary  layer  through  the  Mach  number 
dependent  normal  shock  amplification  factor  derived  by  Ribner  (Reference  16).  This  has  a 
perfect  gas  limit  of  1.8  for  hypersonic  flow  diminishing  to  a  value  of  1.6  at  the  point  of  closest 
approach  to  the  wall  reported  by  Wallace.  For  purposes  of  the  present  comparison  we  chose 
an  average  amplificatioi  factor  of  1.7.  Thus  for  hypersonic  boundary  layers  we  might  expect 
that  for  pure  vorticity  fluctuations  ahead  of  the  probe  shock  a  relationship  of  the  form: 

would  properly  relate  the  Sutera  parameter  A  to  the  turbulence  intensity. 
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In  the  absence  of  detailed  measurements  of  |  v|/V  for  ahypersonic  nozzle  boundary  layer, 

we  chose  to  consider  simply  the  Klebanoff  (e.g.,  Ref.  29*)  values  of  longitudinal  turbulence 
intensity.  From  Figure  2  of  the  cited  reference  we  find  that  a  value  of  |  v  |/V  =  0.15/ 

1.7  =  0.088  occurs  at  y/  8  =  0.021  (also  at  y/8  =  0.002  but  no  heat  transfer  measurements 
were  recorded  in  the  shock  tunnel  sublayer).  Now  from  Figure  6  we  find  for  the  shock  tunnel 

operating  conditions  (M<x>  =  7.9,  T  =  1800#R,  R  /ft  =  2.0  —  10.0  x  106)  that  E  =  0.06  and 

for  A  =  3.8,  q/q^  =  2.25.  This  is  of  the  order  of  the  measured  values  of  excess  heating  at 

y/  8  =  0.05  for  a  number  of  shock  tunnel  runs. 

Finally  a  curve  of  A  vs  y/8  obtained  from  Equation  44  and  the  Klebanoff  profile  as  well  as 
a  curve  of  vs  y/8  obtained  from  the  curve  of  A  vs  y/8  and  Figure  6  appear  as  Figure  7. 

Included  are  the  values  of  q/q^  obtained  from  the  shock  tunnel  measurements.  It  should  be 

noted  that  considerable  scatter  of  experimental  data  prevents  a  completely  satisfactory 
comparison.  However,  one  sees  that  the  unexpectedly  high  heating  rates  measured  (in  some 
cases  over  100%  increase)  and  their  diminution  with  increased  y/8  are  on  the  whole  predict¬ 
able.  K  may  therefore  be  concluded  that  this  analysis  of  stagnation  zone  heating  is  capable  of 
predicting  the  high  heating  rates  encountered  when  the  free  stream  is  turbulent. 

2.  ARC-HEATED  WIND  TUNNEL  TURBULENCE 

An  in-house  program  to  characterize  the  turbulence  in  a  typical  air  arc-heated  hypersonic 
wind  tunnel  is  currently  under  way  at  the  Air  Force  Flight  Dynamics  Laboratory.  The  initial 
phase  of  the  program  involves  measurement  of  as  many  forms  of  unsteady  behavior  in  the 
tunnel  as  possible  with  existing  instrumentation  and  ultimate  auto  and  cross  correlation  of  the 
data.  The  analysis  has  not  been  completed;  however,  certain  of  the  preliminary  results 
illustrate  the  basic  nature  of  the  turbulence  problem  discussed  in  Section  II. 

Basically  the  experiments  involved  introducing  a  Kaman  model  K1808  acoustic  transducer, 
with  nearly  flat  frequency  response  from  10  Hz  to  20,000  Hz  and  pressure  range  from  120  to 
200  db  sound  pressure  level  (SPL),  mounted  in  either  a  cylindrical  or  rectangular  flat-faced 
probe  into  the  test  effluent  of  the  AFFDL  4-megawatt  electrogaadynamics  facility  (EGF) 
(Reference  30).  (Also  see  Figure  1.)  In  order  to  prevent  blocking  by  the  1  inch  x  6  inch 
rectangular  flat-faced  stagnation  probe,  an  18-inch  diameter  nozzle  was  utilized.  Measure¬ 
ments  of  heater  voltage,  pressure,  microphone  output,  as  well  as  the  output  from  an  RCA  6217 
photomultiplier  looking  directly  at  the  gas  cap  radiation  were  recorded  via  an  Ampex  CP  100 
FM  tape  recorder  with  a  frequency  range  of  0  to  27  kHz, 

Of  particular  interest  are  the  microphone  and  photomultiplier  outputs.  A  sample  of  the  tape 
records  is  presented  in  Figure  8  and  the  corresponding  power  spectra  follow  in  Figures  9a 
and  9b.  One  sees  that  a  rather  broad  spectrum  of  turbulence  exists  at  the  axisymmetric 
probe,  stagnation  point.**  E  is  observed  from  the  scope  traces  that  the  microphone  and 
photomultiplier  display  some  degree  of  cross  correlation  as  would  be  expected.  This  follows 
from  the  fact  that  under  the  conditions  of  the  test  the  upper  frequency  limit  of  the  microphone 
(20  kHz)  corresponds  to  a  wavelength  in  the  gas  cap  region  of  about  0.2  ft.  Thus  over  the 


*  In  Reference  29,  Figure  2,  the  component  of  turbulence  intensity  in  the  flow  direct! 'a  is 
labeled  u’/Uo  and  in  our  notation  this  is  |v|/VQ. 

**  Very  similar  power  spectral  densities  occur  for  the  two-dimensional  probe  tests. 
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Figure  7.  Variation  of  Stagnation  Point  Heat  Transfer  Through  a  Turbulent  Boundary  Layer 
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Figure  8.  Photomultiplier  and  Acoustics  Probe  Oscilloscope  Traces 
(Experimental  Configuration  Corresponds  to  Figure  1) 

range  of  frequencies  examined,  the  data  corresponds  to  gross  property  variations  within  the 
entire  shock  layer.  The  same  considerations  based  on  the  foregoing  analysis  place  the  cutoff 
frequency  fQ  at  about  200  kHz.  This  is  based  on  Sutera’s  cutoff  wavelength  approximately 

equal  to  2.6  times  the  classical  Hlemenz  boundary  layer  thickness  (Reference  25)  which  in 
turn  is  set  equal  to  the  estimated  probe  face  boundary  layer  thickness  of  0.1  inch.  It  can  thus 
be  stated  that  velocity  disturbances  normal  to  a  stagnation  line  with  frequencies  below  200  kHz 
will  be  amplified  upon  passing  through  the  boundary  layer.  The  dissipation  mechanism,  of 
course,  will  predominate  at  higher  frequencies. 


3.  STAGNATION  HEATING  IN  PLASMA  JET 

Recent  measurements  of  stagnation  point  heating  have  been  made  is  a  high  enthalpy 
turbulent  nitrogen  jet  exhausting  to  atmosphere  (Reference  31).  The  jet  was  produced  by  a 
l.5-megawatt  arc  heater  and  was  caused  to  impinge  normally  on  a  large  instrumented  flat 
plate  positioned  at  several  axial  stations.  Nominal  operating  conditions  reported  were  :  jet 
exit  radius  r  =  0.375  inch,  exit  Mach  number  *  0.782,  and  exit  stagnation  enthalpy  *  4460 

BTU/lb  giving  an  exit  stagnation  temperature  of  5000*K.  The  jet  How  has  been  previously 
investigated  and  all  important  thermodynamic  properties  are  reported  in  Reference  32. 

The  measured  stagnation  point  heat  transfer  rates  reported  in  Reference  31  were  compared 
with  the  values  expected  based  on  the  theory  of  Fay  and  Riddell.  The  measured  pressure 
distribution  on  the  plate  was  employed  to  establish  the  required  stagnation  point  velocity 
gradient.  This  comparison  revealed  that  the  measured  stagnation  heat  flux  was  between  two 
and  three  times  greater  than  the  laminar  theory  predicts  and  it  was  suggested  that  this  dis¬ 
crepancy  was  due  to  the  presence  of  turbulence  in  the  jet. 


Previous  investigations  of  the  jet  via  high  speed  motion  pictures  revealed  relatively  little 
instability.  Measurement  of  sound  pressure  level,  however,  revealed  a  120-db  white  noise 
spectrum  over  a  range  from  20  Hz  to  10,000  Hz.  This  result  may  be  compared  with  the  present 
findings  for  the  hypersonic  jet  flow  case  (see  ARC-HEATED  WIND  TUNNEL  TURBULENCE 
in  this  Section)  where  again  a  white  noise  spectrum  is  observed. 
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Figure  9a.  Power  Spectral  Density  Analysis,  Acoustic  Probe 
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Figure  9b.  Power  Spectral  Density  Analysis,  Photomultiplier 
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In  an  attempt  to  assess  quantitatively  the  possible  effects  of  free-stream  turbulence  on  the 
increased  stagnation  heating,  the  authors  assumed  (in  the  absence  of  actual  turbulence  level 
measurements)  that  the  axial  distribution  of  turbulence  intensity  was  similar  to  that  found  by 
Corrsin  and  Uberol  (Reference  33).  S  was  assumed  that  for  a  particular  value  of  me.  n 
velocity,  nondime  naionalized  with  respect  to  the  exit  velocity,  the  turbulence  level  would  be 
that  predicted  by  Corrsin  and  Uberoi  for  the  same  nondimensional  mean  velocity.  Since  the 
axial  mean  velocity  in  plasma  Jet  decayed  more  rapidly  than  typical  cold  jet  results,  the 
foregoing  assumption  leads  to  a  somewhat  higher  turbulence  level  at  a  given  axial  location  in 
the  jet. 

Both  Corrsin  and  Uberoi’s  results  for  the  cold  jet  (u/ue,  vs  x/re)*  and  corresponding 

results  for  the  plasma  jet  are  plotted  in  Figure  10.  Condition  (1)  in  the  figure  corresponds  to 
the  plasma  jet  results  whereas  Condition  (2)  refers  to  Corrsin  and  Uberoi’s  cold  jet  measure¬ 
ments. 

These  results  may  be  utilized  in  conjunction  with  the  present  theory  to  predict  the  turbulent 
stagnation  heating  rates.  The  procedure  is  the  same  as  that  outlined  in  this  Section  under 
SHOCK  TUNNEL  TURBULENT  NOZZLE  BOUNDARY  LAYER.  Equation  44 relates  the  param¬ 
eter  A  to  the  measured  turbulence  level.  The  Eckert  number  was  again  obtained  from  Equa¬ 
tion  42  and  was  found  to  vary  from  0.12  at  the  exit  to  0.03  at  the  last  data  station.  Values  of 
q/q^  were  calculated,  employing  Figure  6,  corresponding  to  the  turbulence  levels  given  by 

Conditions  (1)  and  (2).  This  produced  two  heating  curves  which  are  plotted  along  with  the 
experimental  results  in  Figure  10.  It  appears  that  the  experimental  results  are  well  bracketed 
by  the  theory  \f>  to  about  25  nczsle  exit  radii  downstream.  Beyond  this  point  it  might  well  be 
expected  that  errors  in  heat  transfer  measurement  and  plate  pressure  distribution  would 
accumulate  rapidly. 


*  T  is  the  turbulence  intensity,  u andr  are  the  nozzle  exit  velocity  and  radius,  respectively. 

u  O  O 
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Figure  10.  Stagnation  Point  Heat  Transfer  Rate  Along  the  Centerline 
of  a  Nitrogen  Plasma  Jet 


AFFDL-TR-67-195 


SECTION  V 

CONCLUDING  REMARKS 


The  analytical  approach  discussed  in  this  report  appears  to  provide  a  means  of  assessing 
quantitatively  the  effects  of  velocity  disturbances  ahead  of  the  stagnation  point  boundary  layer 
on  heat  transfer.  There  are  several  other  boundary  layer-disturbance  interaction  problems 
that  must  be  explored  along  with  a  broad  class  of  shock-turbulence  Interaction  problems  in 
order  to  fully  understand  the  role  of  fre e-stream  disturbances  on  hypersonic  stagnation  zone 
heat  transfer.  As  a  con*  riuence  of  the  favorable  agreement  between  the  present  analysis  and 
the  shock  tunnel  boundary  layer  experiments,  a  number  of  tentative  conclusions  can  be  drawn: 

1.  The  turbulent  intensity  profile  in  the  hypersonic  boundary  layer  investigated  may  be 
approximated  by  a  Klebanof  distribution.  The  maximum  turbulence  intensity  in  the  boundary 
layer  would  then  be  of  the  or  .er  of  10%. 

2.  Amplification  at  longitudinal  turbulence  intensities  via  a  normal  shock  Interaction  process 
may  be  given  in  the  high  Mach  number  limit  by  Ribner’s  shock-shear  wave  interaction  analysis. 

3.  A  stagnation  point  heat  transfer  probe  with  dimensions  sufficiently  small  with  respect  to 
the  smallest  energy  bearing  eddys  might  be  employed  successfully  to  measure  free-stream 
or  boundary  layer  turbulence  levels. 

The  arc-heated  wind  tunnel  results  are  mostly  qualitative  but  reveal  most  conclusively  that 
the  effluent  contains  turbulent  energy  as  evidenced  by  the  white  noise  spectrum.  No  quanti¬ 
tative  estimates  of  turbulence  intensity  can  be  offered  at  this  time  and  must  remain  the  purpose 
of  further  exploration  both  in  the  core  and  shroud  flow  regions.  An  order  of  magnitude  analysis 
involving  typical  eddy  sizes  suggests  that  velocity  disturbances  normal  to  a  stagnation  line 
with  frequencies  up  to  200  kHz  will  be  amplified  upon  passing  through  the  boundary  layer. 
The  dissipation  mechanism  will  predominate  at  higher  frequencies. 

A  possibility  exists  to  create  stagnation  heating  rates  of  the  order  of  twice  the  laminar 
value  in  an  arc-heated  wind  tunnel  simulation  of  high  heat  transfer-high  shear.  Such  a  condition 
corresponds  to  condition  “C”  in  Figure  6  where  the  Eckert  number  is  about  the  same  as  for 
the  reported  shock  tunnel  experiments.  This  might  be  accomplished  by  rendering  the  effluent 
turbulent  in  some  manner  making  sure  that  the  longitudinal  velocity  turbulence  level  is  of  the 
order  of  10%. 

The  plasma  Jet  experiments  have  served  to  further  corroborate  the  analytical  results.  Here 
again,  the  actual  turbulence  level  was  not  measured,  but  a  reasonable  assumption  is  that  it 
corresponds  to  measured  cold  jet  results. 
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APPENDIX 

LISTING  OF  SPECTRAL  INTERACTION  FUNCTIONS 


In  the  following  expressions  i  is  a  dummy  index  and  quantities  with  subscripts  n-i  or  i-n 
are  to  be  replaced  with  zero  whenever  their  subscripts  are  negative  or  zero. 


Corresponding  to  Equation  12  we  list  the  I 
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Corresponding  to  Equation  35  we  obtain  the  following  I 


X(  i.i) 

=  1,  =0 
(  2  ,l) 

I3  0) 

=  2ul(2)U2  (2) 

lA.  . 

*\  1 1 

=  2  Mai  V  2  {  2) 

1 » (i ) 

=  2  Wl  (2)  W2  (2) 

29 


AFFDL-TR-67-195 


*  wlu),.(.)  +  <(«)  vui  +w'2(2)v,(2)  ~‘rwM*)V> 


’  (  7,  I  ) 


=  k|  (^1(2)  V2(2)"*”Vl(a)  V2  (2)+  4  V2(2)VI  (2)  ~  2  V2  (2)  V  I  (2)  * 

=  J_  /  Wife)  *2(l)  +  S(2)  W2(2)  X 

l  *  '  r>  o  ' 


The  spectral  Interaction  functions  I  0.  spearing  in  Equation  36  are: 
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Finally  the  I  ^  jj,  I  ^  ^  corresponding  to  Equation  38  and  the  I  ^  ^  2y  I  (2,  2)  corresPond“ 
ing  to  Equation  39  are  listed. 

I(i,i)  =  v2{2)  ^i(i)  +  vi(2)fl2  (i) 
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